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Abstract 

In  this  contribution,  a  dynamical  nonlinear  model  structure  for  the  identification  of  the  electrical  behavior  of  a  solid  oxide  fuel  cell  is  developed 
using  expert  knowledge  from  material  scientists.  The  estimation  of  the  model  parameters  can  be  realized  by  an  iterative  algorithm  based  on 
analytical  equations  to  determine  parameter  estimates  in  each  iteration  step.  An  enhancement  of  the  estimation  can  be  achieved  by  an  additional 
numerical  optimization  of  the  model  parameters.  Furthermore,  this  model  can  be  extended  to  represent  the  time  dependency  of  the  model  due  to 
deterioration  of  the  fuel  cell  materials. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

For  many  purposes  a  mathematical  model  of  a  solid  oxide 
fuel  cell  (SOFC)  is  required.  For  instance,  this  model  allows  to 
separate  the  formation  effects  from  the  usual  current/voltage- 
dynamics  during  the  first  startup  procedure  of  a  SOFC  [1]. 
During  operation,  it  could  be  used  to  develop  a  control  strategy 
for  load  changes.  And  finally,  a  dynamical  model  can  represent 
the  effects  of  deterioration  of  a  SOFC-system.  Recent  publica¬ 
tions  like  [2-4].  address  the  modeling  of  these  systems.  In  [2] 
and  [4]  a  physical  approach  was  used  to  represent  the  dynamic 
behavior  of  the  fuel  cell.  To  obtain  a  lower  order  linear  model 
of  the  SOFC-dynamics  an  identification  procedure  was  pro¬ 
posed  in  [3].  The  paper  is  organized  as  follows:  in  Section  2,  an 
appropriate  model  structure  will  be  developed.  Afterwards,  an 
identification  for  a  time-invariant  parametric  model  is  proposed 
in  Section  3.  In  order  to  improve  the  mathematical  representa¬ 
tion,  a  time-variant  model  will  be  estimated  in  Section  4.  Finally, 
in  Section  5  a  simple  approach  to  represent  degradation  effects 
is  given. 
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2.  Development  of  a  model  structure 

In  this  paragraph  a  model  structure  for  the  current- 
density/voltage-loss-behavior  will  be  developed.  The  cathode 
material  of  the  considered  SOFC  is  (La,Sr)Mn03_  which  is 
abbreviated  as  LSM.  Yttrium  doped  Zirconia  8YSZ  (8mol% 
Y2O3  doped  ZKE)  is  the  electrolyte  material.  Finally,  the  anode 
consists  of  a  Ni/8YSZ-cermet. 

In  Fig.  1,  measurement  data  of  a  SOFC  are  given.  Based  on 
the  same  measurement,  the  corresponding  voltage-loss  versus 
current-density  is  depicted  in  Fig.  2.  The  latter  plot  shows  that 
the  electrical  behavior  is  dominated  by  a  static  transfer  element 
and  that  dynamic  effects  are  less  significant.  This  suggests  the 
following  modeling  approach 

u(kT)  =  ux(kT)  +  un(kT),  (1) 

uN(kT)  =  N(jc(kT)),  (2) 

uH(kT)  =  Z-1  {ife-1)}  *  <P(jc(kT)),  (3) 

where  u(kT)  is  the  voltage  loss  and  jc(kT)  denominates  the  cur¬ 
rent  density  of  the  fuel  cell,  T  is  the  constant  sampling  time  and  k 
is  the  index  for  the  sampled  measurement  data.  In  ( 1 ),  the  voltage 
loss  is  the  sum  of  a  nonlinear  function  N(jc)  of  the  current  density 
jc  and  the  response  of  a  Hammerstein-model  [5]  consisting  of  a 
linear  system  L(z_1)  whose  input  signal  is  a  function  0(/c)  of  the 
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Fig.  1.  Current-density/voltage-loss-measurements  vs.  time. 


|  [Cathode  grain  |  [Triple  phase  boundary  |  [Electrolyte 


Fig.  3.  TPB  in  equilibrium  before  a  current  density  change. 


Fig.  4.  TPB  immediately  after  an  increase  of  the  current  density. 


current  density  jc .  The  symbol  in  (3)  indicates  the  convolu¬ 

tion  operator.  The  current-density/voltage-loss-measurement  in 
Fig.  2  is  mainly  represented  by  the  function  N(jc).  Some  assump¬ 
tions  have  to  be  made  in  order  to  specify  the  model  structure 
(l)-(3)  in  more  detail.  It  is  known  that  two  reversible  processes 
in  the  materials  are  responsible  for  the  dynamic  behavior  of  the 
SOFC  which  will  be  represented  by  the  linear  system  L(z~ 1 )  and 
the  unknown  function  0(/c): 


triple  phase  boundary  is  in  the  equilibrium  state  for  the  new 
higher  current  density  (Fig.  5).  Thus,  the  conductivity  is 
increased  by  the  increase  of  the  triple  phase  boundary. 

The  model  order  of  the  linear  system  L(z~l )  is  chosen  as 
ul  =  2,  in  order  to  represent  the  two  considered  dynamical 
effects,  leading  to  the  transfer  function 


( 1 )  Temperature  effect:  a  changing  current  density  causes  a  tem¬ 
perature  change  of  the  material  which  results  in  a  decrease 
or  an  increase  of  conductivity.  The  time  constant  of  this 
process  is  about  3-5  min. 

(2)  Structural  effect:  a  changing  current  density  causes  a 
change  of  the  local  oxygen  partial  pressure  at  the  cath¬ 
ode/electrolyte  interface.  The  oxygen  partial  pressure  has 
an  impact  on  the  oxygen  vacancy  concentration  in  the  cath¬ 
ode  grain  [6].  Therefore,  the  size  of  the  electrochemical 
active  triple  phase  boundary  at  the  cathode/eletrolyte  inter¬ 
face  is  affected.  In  Fig.  3,  the  triple  phase  boundary  at 
the  cathode-electrolyte  interface  in  equilibrium  is  depicted 
schematically.  An  instantaneous  current  density  step  causes 
an  increase  of  the  area  specific  current  density  at  the  triple 
phase  boundary.  This  effect  is  shown  in  Fig.  4.  The  size  of 
the  triple  phase  boundary  increases  due  to  the  higher  area 
specific  current  density.  This  effect  could  be  explained  by 
a  depletion  of  oxygen  in  the  cathode  grain.  After  1-5  h,  the 


,,  _k  bo  +  biz  1  +b2Z  2 

L\Z  )  —1  i  —2  ' 

1  +  a\z  1  +  ci2Z  z 

In  order  to  separate  the  stationary  behavior  and  the  dynami¬ 
cal  behavior,  the  linear  systems  output  signal  has  to  vanish  for 
constant  input  signals.  For  constant  input  signals  jc,  the  output 
signal  is  solely  determined  by  the  nonlinear  function  N(jc).  This 
leads  to  a  constraint  for  the  numerator  of  the  linear  system.  To 
meet  this  condition,  the  linear  system  L(z~l)  has  to  have  a  dif¬ 
ferentiating  behavior.  Time  discrete  systems  with  a  zero  at  z  =  1 
show  this  property.  The  constraint 

bo  +  b\  +  Z?2  =  0  (4) 

for  the  model  parameters  bo,  b\  and  Z?2  guarantees  that  the  linear 
systems  response  equals  zero  in  steady- state. 

It  is  also  known  that  the  time  constants  of  the  two  processes 
depend  on  the  sign  of  djc/dt  =  jc.  Therefore,  the  estimation 
of  two  different  linear  systems  is  necessary,  one  system  for 
increasing  current  densities  and  a  different  one  for  decreasing 
current  densities.  For  jc  >  0  the  linear  model  L(z~l)  is  denoted 
as  Lup(z-1)  and  for  jc  <  0  it  is  denoted  as  Tdownfe-1)- 


Fig.  2.  Current-density  vs.  voltage-loss. 


Fig.  5.  TPB  in  equilibrium  after  the  current  increase. 
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Fig.  6.  Structure  of  the  identification  model  in  most  cases. 

For  the  nonlinearity  N(jc),  a  polynomial  approach 

NUc)  =  c0  +  Cijc  +  C2j'c  +  •  •  •  +  cnNj”N  (5) 

is  proposed.  The  degree  n^  =  3  for  the  polynomial  N(jc)  in  (5)  is 
adequate  in  most  cases. 

This  model  structure  is  depicted  in  Fig.  6.  In  the  sequel,  the 
unknown  model  parameters  of  this  model  structure  will  be  deter¬ 
mined  from  measurement  data. 

3.  Iterative  identification  of  the  model  parameters  of  the 
time-invariant  model 

In  this  section,  a  parameter  identification  algorithm  for  the 
time-invariant  model  structure  in  Fig.  6  is  presented.  By  analysis 
of  measurement  results,  it  turns  out  that  the  unknown  func¬ 
tion  0(/c)  could  be  chosen  as  <p(jc)  =  ft-  Therefore,  the  model 
parameters  of  the  time-invariant  model  structure  in  Fig.  6  are 


$TI  —  I/O?  Oi,  C 2,  C3;  <21, up,  <72, up?  <7l,down?  <72, down? 

^0,up?  ^l,up?  ^2, up?  ^0,down?  ^l,down?  ^2, down]-  (6) 

The  indices  “up”  and  “down”  of  some  model  parameters  in 
(6)  denote  their  affiliation  to  the  linear  submodels  Lup  (z-1)  or 
^down  (z  !)• 

The  basic  approach  in  parameter  identification  is  to  formu¬ 
late  the  identification  problem  as  an  optimization  problem.  Most 
commonly,  a  quadratic  performance  index 

N 

•^Tl(^Tl)  =  ^  ^(jmeasC^)  —  JmodeK^TI’  ^))  (7) 

k=  1 

is  defined,  where  ymeas(£)  are  measured,  noisy  data  from  the  out¬ 
put  of  the  process  and  ymodei(^Ti?  k)  are  the  simulated  output  data 
from  the  used  identification  model  (Fig.  6)  with  the  parameter 
vector  0j\.  A  numerical  optimization  algorithm  could  minimize 
the  function  (7)  under  the  constraint  (4).  With  sufficient  initial 
values  for  this  optimization,  the  determination  of  the  global  min¬ 
ima  of  the  function  /ti(^ti)  is  more  likely.  The  determination  of 
these  initial  values  is  given  by  an  iterative  algorithm  using  the 
dominating  influence  of  the  static  element  N(jc)  on  the  transfer 
behavior.  The  steps  of  this  algorithm  are  depicted  in  the  function 
chart  in  Fig.  7: 


Fig.  7.  Function  chart  of  the  identification  algorithm. 

(1)  The  difference  signal  w^ident  =  ^meas (kT)  —  uaikT)  is  calcu¬ 
lated  using  the  up  to  date  Hammerstein-Model  unikT).  In  the 
first  initialization,  this  signal  is  set  to  unikT)  =  0,  to  neglect 
the  influence  of  the  small  signal  uuikT)  on  the  absolute  value 
u(k)  (see  Fig.  2). 

(2)  The  model  parameters  6 n  of  the  polynomial  N(jc)  are  esti¬ 
mated  by  the  least-squares-algorithm  [5,7]  using  un, ident 
which  was  computed  in  step  1 . 

(3)  With  the  model  parameters  0^  it  is  possible  to  separate 
the  signals  unikT)  and  u^ikT)  in  the  measurement  data 
Mmeas (kT).  The  signal  uH-l(ient(kT)  =  umeas(kT)  -  uN(kT)  is 
computed. 

(4)  Using  w//;ident(£T)>  the  linear  subsystems  Lup(z_1)  and 
^downfe-1)  will  be  estimated  taking  constraint  (4)  into 
account.  This  estimation  is  computed  with  Eq.  (A.5), 
derived  in  Appendix  A. 

(5)  Using  all  currently  estimated  parameters  gives  the  modeling 
error  of  the  model  parameter  set  J(0ji)  =  Ec2(k). 

(6)  If  J(0ji)  <  /max  is  not  fulfilled,  the  algorithm  jumps  to  step 
1  using  the  estimated  signal  unikT). 

If  <  /max?  a  numerical  minimization  of  J(6ji)  with 

respect  to  the  model  parameters  0j  1  gives  the  improved  esti- 
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Fig.  8.  Voltage  loss  of  the  model  and  measurement  data. 

mated  value  for  0Tl.  This  minimization  is  necessary  because  it 
is  not  guaranteed  that  the  iterative  execution  of  steps  1-5  yields 
the  global  minimum  of 

The  result  of  an  estimation  using  the  measurement  data  from 
Fig.  1  is  shown  in  Fig.  8.  This  plot  shows  that  the  model  struc¬ 
ture  (1)— (3)  is  an  appropriate  modeling  approach,  because  the 
modeling  error  is  very  small.  In  Appendix  B,  a  transformation  is 
given  to  separate  the  two  dynamical  processes.  With  (B.7),  the 
time  constants  of  the  two  processes  can  be  determined,  which 
makes  the  allocation  of  the  two  processes  possible.  For  this  iden¬ 
tification  result,  the  time  constants  are: 

Ti,Up  =  3.52  min  ri,doWn  =  4.08  min 

(o) 

^2,up  =  2.47  h  T2,down  =  4.9  h 

The  values  for  the  time  constants  in  (8)  indicate  that  process  1 
is  the  temperature  effect  and  process  2  is  the  structural  effect.  In 
Fig.  9,  the  voltage  losses  of  these  two  processes  are  shown  for  a 
section  of  the  measurement  period.  The  term  uje  designates  the 
temperature  effect  voltage  and  use  designates  the  voltage  of  the 
structural  effect.  The  plot  shows  that  the  voltage-loss  of  the  struc- 


t/h 

Fig.  9.  Separation  of  the  two  dynamical  effects. 


Fig.  10.  Identified  current-density/voltage-loss-characteristic  N(jc). 


tural  effect  is  by  far  more  significant  for  the  dynamical  behavior. 
An  additional  advantage  of  the  identification  procedure  pro¬ 
posed  in  this  paper  is  that  the  current-density /voltage-loss  char¬ 
acteristic  could  be  estimated  simultaneously  to  the  estimation 
of  the  SOFC-dynamics.  The  estimated  current-density/voltage- 
loss  characteristic  is  shown  in  Fig.  10. 

The  proposed  identification  method  was  also  applied  suc¬ 
cessfully  on  other  measurement  data  taken  from  solid  oxide  fuel 
cells,  resulting  to  different  values  for  the  model  parameters  than 
those  determined  in  this  paper. 

4.  Time-variant  model  identification 

Until  now,  a  time-invariant  parametric  model  was  identified. 
However,  it  is  known  that  the  electrical  behavior  of  SOFC- 
sy stems  changes  during  large  time-periods.  Fig.  1 1  demonstrates 
the  time  variance  of  the  considered  fuel  cell.  The  continuous  line 
is  the  current-density/voltage-loss  measurement  which  was  used 
for  the  identification  of  the  voltage  loss  model  in  Section  3.  The 
measurement  data  between  410  and  490  h  is  depicted  as  a  dotted 
line  in  the  same  figure.  Clearly,  it  can  be  recognized  by  analyz¬ 
ing  this  plot  that  the  current-density/voltage-loss  characteristic 
could  not  be  assumed  as  time-invariant.  The  linear  component  of 
the  characteristic  c\  has  changed  most  significant  during  oper¬ 
ation  of  the  fuel  cell.  Therefore,  a  time-variant  model  for  the 
voltage  loss  will  be  estimated  in  this  section  in  order  to  improve 
the  representation  of  the  real  physical  behavior. 


Fig.  11.  Measurements  of  current-density  vs.  voltage-loss. 
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A  simple  but  common  method  used  in  the  field  of  iden¬ 
tification  is  the  minimization  of  an  exponentially  weighted 
sum 

i 

JTV(0TY,  i)  =  «*(*)  -  Vmodel^TvO'),  V?  (9) 

wp(k) 

of  the  modeling  error.  The  weighting  factor  has  to  weight  the 
present  output  error  stronger  than  the  errors  from  the  past  [8]. 
The  forgetting  factor  >  0  has  to  be  chosen  suitable  to  detect 
variations  of  the  model  parameters.  On  the  other  hand,  the 
forgetting  factor  should  not  be  too  large  in  order  to  keep  the 
influence  of  the  measurement  noise  on  the  identification  result 
limited.  For  each  index  value  i  a  numerical  minimization  of  the 
weighted  performance  index  (9)  gives  a  snapshot  of  the  sys¬ 
tem  behavior  at  the  time  IT,  because  the  vector  Gji(i)  of  the 
model  parameters  is  also  time- variant.  The  numerical  minimiza¬ 
tion  uses  the  previous  estimated  time- variant  parameter- vector 
0Tl(i  —  1)  as  the  initial  value  for  the  estimation  of  the  current 
estimation  of  0TI(/),  because  the  changes  between  two  time 
instants  can  be  assumed  as  very  small.  The  first  identifica¬ 
tion  uses  the  time-invariant  parameters  0Tl  as  initial  values. In 
advance  of  the  execution  of  the  time-variant  identification,  it 
is  necessary  to  decide  which  model  parameters  will  be  time 
variant.  A  high  degree  of  freedom  causes  a  very  good  identifi¬ 
cation  accuracy,  but  it  will  not  be  possible  to  interpret  the  result. 
Therefore,  the  degree  of  freedom  has  to  be  reduced.  As  already 
mentioned,  Fig.  11  demonstrates  that  the  parameter  c\,  which 
represents  the  ohmic  component  of  the  characteristic,  changes 
most  significantly.  Therefore,  it  will  be  sufficient  to  assume  that 
only  c\  varies  with  respect  to  time.  All  other  model  parameters 
will  be  assumed  as  time-invariant.  The  result  of  the  numerical 
optimization  for  the  identification  of  the  time-variant  c\{i)  is 
shown  in  Fig.  12.  The  identification  algorithm  for  the  param¬ 
eter  0Tl  starts  at  i  =  5690.  Fig.  12  shows  that  the  parameter  c\ 
increases  monotonically,  if  the  voltage  loss  is  increasing.  A  con¬ 
stant  or  a  decreasing  voltage-loss  has  a  very  faint  influence  on 


the  change  of  c\.  This  result  is  surprising,  because  it  shows 
that  an  increase  of  voltage-loss  has  a  more  significant  influ¬ 
ence  on  the  degradation  of  the  SOFC  than  constant  high  current 
densities. 

This  surprising  new  result  can  be  explained  by  the  structural 
effect.  In  Figs.  3-5,  this  effect  is  sketched  schematically  for  a  cur¬ 
rent  density  increase.  It  takes  1-5  h  to  get  in  the  equilibrium  state. 
During  this  time  period  the  triple-phase-boundary  increases  and 
causes  a  decrease  of  the  voltage-loss.  During  the  settling  time, 
the  voltage  loss  and  the  area  specific  current  density  at  the  triple- 
phase-boundary  is  very  high.  Area  specific  current  densities  have 
a  strong  influence  on  the  materials  at  the  triple  phase  boundary. 
In  [1],  the  area  specific  current  density  causes  a  decrease  of 
losses  during  the  formation  of  the  cathode/electrolyte-interface 
of  SOFCs.  Here,  the  locally  high  loss  during  the  settling  time 
after  an  increase  of  current  density  causes  degradation  of  the 
SOFC  (increase  of  the  ohmic  losses),  the  efficiency  of  the  cell 
diminishes.  This  observed  degradation  process  will  be  modeled 
in  the  next  section. 


5.  Model  for  the  degradation  effects 

In  order  to  predict  the  status  of  the  fuel  cell  materials,  the 
observed  degradation  effect  has  to  be  represented  by  a  math¬ 
ematical  model.  As  demonstrated  in  Fig.  12,  an  increase  of 
the  current  density  causes  an  increase  of  the  model  parame¬ 
ter  c i,  which  represents  the  ohmic  resistance  in  the  current- 
density/voltage-loss-characteristic  N(jc).  A  possible  explanation 
for  this  increase  is  a  locally  concentrated  high  current  density  at 
the  triple-phase-boundary  at  the  cathode/electrolyte-interface. 
A  possible  mathematical  approach  is  the  use  of  the  voltage  loss 
of  the  structural  effect  use  as  an  input  signal  of  a  discrete-time 
integrator.  To  take  into  account  that  c\  does  not  increase  for 
decreasing  current  densities,  only  positive  values  of  use  should 
have  an  impact  on  the  degradation  in  the  model.  The  result  is  the 
model 


c\{k)  =  c\(k  -  1)  +  mu+E(kT)) 

=  c\(k  —  1)  +  T[do  +  d\u^E(kT))  +  d2(u^E(kT))2 
+  d3(u+E(kT ))3  +  d4(u+E(kT))4\,  (10) 


which  is  used  to  represent  the  degradation  (the  increase  of  c\), 
with 


*SE 


use  for  use  >  0 
0  for  use  <  0 


as  the  input  signal.  T  8(u$E(kT))  is  the  difference 
c\(k)  —  c\(k—  1),  hence  the  8  value  could  be  interpreted 
as  the  velocity  of  the  degradation.  The  determination  of  the 
model  parameters 


Fig.  12.  Identification  of  the  time- variant  parameter  c\(k). 


£ieg  =  Ido  d\  d2  d3  d4\ 
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of  the  degradation  model  (10)  is  done  by  a  numerical  minimiza¬ 
tion  of  the  performance  index 

N 

^deg(^deg)  =  ^  ^(jmeas(^)  —  ymodel($deg’  ^)) 
k=  1 

for  the  output  error  ymeas(£)  —  ymodei  of  the  voltage  loss  model, 
where  c\(k)  is  time- variant  and  given  by  (10).  The  other 
model  parameters  are  the  time-invariant  estimates  from  Sec¬ 
tion  3.  In  Fig.  13,  the  estimated  function  8(u^E)  is  shown 
for  the  measurement  data.  The  interception  point  of  the  plot 
with  the  5-axis  gives  the  increase  of  c\  for  the  steady  state 
operation  (uh  =  u$E  =  0).  For  0  <  u$E  <  1  mV  the  degrada¬ 
tion  velocity  8(u^E(kT))  is  increasing  monotonically.  For  high 
values  u$E  <  lmV,  8  is  decreasing.  In  [1],  it  was  demon¬ 
strated  how  high  area  specific  currents  could  improve  the  cell 
performance  during  the  formation  process.  Thus,  a  possible 
physical  explanation  for  the  smaller  degradation  velocity  8  is 
that  the  formation  of  the  cathode/electrolyte-interface  is  con¬ 
tinued,  due  to  very  high  area  specific  current  densities  for 
<c  1  mV. 


systems  in  order  to  allow  a  long  durability,  because  the  expected 
degradation  of  the  fuel  cell  materials  can  be  predicted  by  this 
model. 
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Appendix  A.  Least-squares-estimation  with  a  linear 
constraint 


Given  is  the  performance  index  function 


j(9)  =  (y-  MUfiy  -  Ml) 

(A.l) 

and  the  linear  equation 

N0_  =  b. 

(A.2) 

Objective  is  finding  a  0  which  minimizes  the  function  (A.l)  and 
meets  the  constraint  (A.2).  This  problem  will  be  solved  by  using 
the  method  of  the  Lagrangian  multipliers.  The  new  function 

L(0)  =  J(0)  +  XT(N0  -  b) 

is  defined,  whose  derivation  with  respect  to  the  variable  9  and 
the  Lagrangian  multiplier  X  gives  the  necessary  conditions 

=  -2  MTy  +  2MtM0  +  NTX  =  0,  (A.3) 

0=0 


dL(0) 

do 


9L(0) 

dX 


=  NO-  b  =  0 


(A.4) 


6.  Conclusions  and  outlook 


to  solve  this  optimization  problem.  Using  the  conditions  (A.3) 
and  (A.4)  gives  the  solution 


A  useful  model  structure  for  the  current-density/voltage- 
loss  dynamic  for  a  single  SOFC  is  developed.  The  model 
structure  of  the  generally  nonlinear  model  consists  of  a  static 
element  for  the  steady- state  equilibrium  characteristic  and  a 
switching  Hammerstein-Modell  for  transitions  between  differ¬ 
ent  equilibrium  states.  A  simple  and  fast  identification  algo¬ 
rithm  is  developed  to  identify  the  model  parameters  of  this 
time  invariant  model.  However,  for  long  term  simulations  of 
the  voltage-loss  of  a  fuel  cell,  the  identified  time  invariant 
model  is  not  suitable.  This  can  be  explained  by  changes  in 
the  cathode/electrolyte-interface  caused  by  a  locally  high  area 
specific  current  density.  In  order  to  represent  the  time  depen¬ 
dency,  a  time- variant  model  is  identified  using  an  exponentially 
weighted  performance  index.  A  numerical  optimization  gives 
the  estimated  value  for  a  significant  time  dependent  parameter. 
Using  the  time  dependency  of  the  voltage-loss  model,  a  sim¬ 
ple  approach  is  given  to  represent  the  coherence  leading  to  the 
fuel  cell  degradation  during  fast  changes  of  the  current  den¬ 
sity.  A  low  order  nonlinear  degradation  model  is  estimated  by 
a  numerical  optimization  algorithm.  This  model  can  be  used 
to  develop  control  strategies  for  a  safe  operation  of  SOFC- 


0  =  PM_Ty  +  PN_TQ(b  -  NPMTy)  (A.5) 

if  the  matrices 

P  =  ( mJm)~ 1  and  Q  =  (A7’M)_1A7’)_1 
exists. 


Appendix  B.  Separation  of  the  dynamic  effects 


The  transfer  function  of  the  linear  system 
bo  +  b\z~x  +  b2Z~ 2 


L(z  *)  -  _2, 

1  +a\z  -  1  +«2 z 

has  to  be  transformed  on  the  form 


bo  +  b\+b2  =  0  (B.l) 


1 


1 


L(z  l)  =  hf — - — r +fe  — 
1-Piz-1  1  - 


,-i 


P2Z 


-1 


(B.2) 


L2(z~ ‘) 


where  the  two  dynamic  effects  L\{z  1 )  and L2 (z  1 )  can  be  inves¬ 
tigated  separately.  The  representation  of  (B.2)  is  not  unique, 
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because  two  solutions  are  possible  due  to  the  permutability  of 
the  parameters.  Thus,  one  possible  solution  is 


k\  —  bo  + 


h 


b2  ~  bo  ((-ai/2)  +  ^(a\/4 )  - 


-(ai/2)  +  ^(a\/4)-a2  -  (a2)/  ^-(ai/2)  +  yj (a\ /4)  - 


bo  (^-{a\/2)  +  yj  (a  |  /  4)  -  -  ^2 


-(ai/2)  +  yj («[  / 4)  -  a2  -  ^2/  ^-(ai/2)  +  yJ(a\/4)  - 


Pi  = 


(22 


-(ai/2)  +  yj(a\/4r)-a2 


a\ 


The  time  constants  of  the  two  processes  are 
T 


Ti  = 


In  pi 
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